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Abstract 

We analyze the level crossing rate (LCR) and the average fade duration of the output signal- 
to-noise-ratio (SNR) in generalized switched diversity systems. By using a common approach, we 
study these higher order statistics for two different kinds of configurations: (1) Colocated diversity, i.e. 
receiver equipped with multiple antennas, and (2) Distributed diversity, i.e. relaying link with multiple 
single-antenna threshold-based decode-and-forward (DF) relays. In both cases, we consider the switched 
diversity combining strategies Selection Combining and Switch & Stay Combining (SSC). Whenever 
using threshold-based techniques such as DF or SSC, the output SNR is a discontinuous random process 
and hence classic Rice approach to calculate the LCR is not applicable. Thus, we use an alternative 
formulation in terms of the one and two-dimensional cumulative distribution functions of the output 
SNR. Our results are general, and hold for any arbitrary distribution of fading at the different diversity 
branches. Moreover, we develop a general asymptotic framework to calculate these higher order statistics 
in high mean SNR environments which only needs of the univariate probability density function. 
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I. Introduction 

Switched diversity techniques are trending upward since the introduction of distributed coop¬ 
erative diversity [1]—[3]. In this new setup, a source and a destination are willing to communicate 
with the help of N single-antenna relays. Hence, these intermediate relay stations are regarded 
as a distributed agent, and different combining strategies are employed. These strategies are 
inspired on their single-link multi-antenna receiver counterparts 1 , allowing for high-performance 
with low-complexity terminal equipment. Switched techniques arouse an special interest in this 
context due to their capability of achieving full-order spatial diversity [5]—[8], while all-participate 
strategies entailed loss of spectral efficiency as they require for orthogonal transmission. 

Selection combining (SC) consists in switching each time to the best branch available, usually 
in terms of signal-to-noise-ratio (SNR), to perform the communication between source and 
destination [9, Section 9.8]. In the classical SC colocated multi-antenna receiver, selecting the 
branch with the best instant SNR implies that the N available branches need to continuously 
be monitored. For this reason, Switch & Stay Combining (SSC) was introduced as a way to 
avoid this requirement while allowing for an even simpler receiver. In SSC, a given branch is 
selected as long as the SNR on that branch remains above a given threshold [9, Section 9.9]. 
Nevertheless, as technology has developed, monitoring the SNR on every antenna is no more 
an issue and SSC had been left in low esteem until the appearance of the distributed diversity 
concept. 

Since the introduction of distributed cooperative diversity, counterparts of these classic com¬ 
bining strategies can be found in recent literature. Distributed SC is known since its introduction 
in [10] as Opportunistic Relaying (OR). When considering OR, monitoring the SNR on every 
branch is no longer a matter of equipment’s hardware complexity, but of amount feedback 
required to monitor the end-to-end SNR through every branch and to inform each relaying node 
whether it must retransmit or remain idle. As a result of this concern, SSC experiences a renewal 
of interest in distributed diversity schemes, since feedback is only required to request a switch to 
another relaying branch. To the best of our knowledge, the concept of distributed SSC was first 
introduced in [11], [12]. These distributed counterparts have fostered numerous works devoted 

'For the sake of clarity, we will refer to the single-link multi-antenna diversity counterpart as colocated diversity using the 
notation introduced in [4], 
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to characterize the performance as much of OR, [2], [7], [8], [10], [13], as of distributed SSC 
systems [4], [11], [12], [14]—[17]. Despite this evident interest, most analyses are focused on 
first order metrics such as the outage probability. 

Second order statistics such as the level crossing rate (LCR) or the average fade duration (AFD) 
provide useful information related with the rate of change of a stationary random process, thus 
providing additional information about the dynamics of the process. Specifically, the LCR denotes 
how often the random process crosses a given threshold, whereas the AFD states the average 
amount of time that the random process remains below that threshold level. The knowledge of 
these second order statistics finds a variety of applications in the modeling and design of wireless 
communication systems, such as the amount of feedback required to successfully perform the 
communication or the latency of a link. The pioneering work by Rice [18] has allowed the 
computation of these metrics in a very general fashion, in terms of the joint distribution of the 
random process and its time derivative. One of the advantages of Rice’s approach is its flexibility, 
as apparently mild conditions are imposed to the random process: stationarity, continuity, and 
smoothness of the correlation coefficient. 

While these conditions hold for most common fading scenarios and hence the LCR and AFD 
have been already studied for different fading distributions in diversity combining scenarios [13], 
[19]—[22], the assumption of continuity in processes derived from threshold-based techniques, 
such as threshold-based decode-and-forward (DF) or SSC, is not realistic and hence Rice’s 
framework finds an obstacle: the output SNR of dual-hop link with DF relaying and also the 
output SNR of a SSC scheme, either colocated or distributed, are inherently discontinuous. For 
this reason, analyses dealing with the LCR characterization of threshold-based techniques are 
scarce. To the best of our knowledge, the LCR characterization of the SSC scenario has only been 
studied in [23] under zero temporal correlation assumption using Rice’s approach. Its results can 
be regarded as approximated: their accuracy was not corroborated by simulation and its plotted 
curves do not follow the shape of ours, which are corroborated by simulation. Nevertheless, the 
threshold-based DF obstacle has been recently circumvented by overlapping the results obtained 
from Rice’s approach in different continuous crossing events [13]. 

In this paper, we calculate in a general fashion the LCR and the AFD of processes derived 
from switching combining techniques using an alternative approach for calculating the higher 
order statistics of sampled processes [24], which does not impose continuity as a condition. 
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Furthermore, DF and switched diversity combining systems are implemented in a discrete time 
fashion, hence our processes are sampled, thus this other approach is the suitable one for our 
interest. We express the LCR and the AFD in terms of the univariate and bivariate CDF of the 
output SNR. 

We investigate the SC, OR and SSC techniques in arbitrary fading conditions in this sampled 
fashion. Our results are valid both for the colocated and distributed configurations, and hold 
for any arbitrary fading distribution. We regard the independent and identically distributed (IID) 
branches and independent non-identically distributed (InID) case, allowing for arbitrary temporal 
correlation model of the per-branch fading for every scenario except for the colocated and 
distributed InID SSC case, where unfortunately we can only consider zero temporal correlation, 
i.e. independence between samples of the branch fading envelope. Furthermore, through our 
analysis, we find that the SSC technique in time-correlated fading environments benefits from 
having more than two diversity sources, against the usual recommendation which states that it 
does not [25]. 

Finally, but not less interesting, we also derive simple approximated expressions in terms of 
the fading univariate probability density function (PDF) of the LCR and AFD in high mean 
SNR fading environments. That is, we achieve to accurately calculate this higher order statistics 
for SNR levels much lower than its mean value by only using the one-dimensional PDF of the 
process. 

The remainder of this paper is organized as follows. Section II introduces the scenarios of 
interest. In Section III we present the general analysis for the LCR and AFD in the investigated 
scenarios. Then, in Section IV, we discuss the asymptotic behavior of the LCR and AFD for high 
mean SNR. Section V shows numerical results which give validity to our theoretical expressions. 
Finally, the conclusions of our work are discussed in Section VI. 

II. System Model 

We will consider two different spacial diversity configurations, regarded as colocated diversity 
and distributed diversity. Colocated diversity is the conventional multi-antenna receiver [4], where 
the diversity sources are gathered together at the receiver; thus they are fully available all the 
time 2 . Figure 1.A shows a colocated diversity scenario with two branches, N = 2, where Z t [n] 

2 We will discuss the diversity sources availability at the receiver side further on. 
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with i G [0... N — 1] represents the SNR in the nth time interval on the /1h branch, i.e. /1h 
antenna. In the distributed case, for the sake of spectral efficiency, the diversity sources are not 
available all the time at the receiver side, hence the receiver must request them to the relay 
stations. Figure l.B shows a distributed diversity scenario with two branches, N — 2, where 
Z, ] [n], Zj 2 [n] with i G [0 ... N — 1] represent the SNR on the first and second hop of the /1h 
branch, respectively, and Zj [n] with i G [0 ... N — 1] is the overall SNR through the zth branch, 
i.e. ith relayed channel. In both configurations SC or SSC is performed to get a system output 
with SNR Z [n]. 


A. 



Z 0 [n] 

V\ 



Fig. 1. A. Colocated diversity scenario. B. Distributed diversity scenario. 


On each of these scenarios, we assume that additive white gaussian noise (AWGN) affects ev¬ 
ery branch, and that the stationary processes that characterize the different links are independent. 
We also consider as in [13] that switching from one branch to another is to be performed every 
time interval Tg, short enough so that the instantaneous SNR on each antenna remains constant 
through Tg, nevertheless it will change from an interval to the next one. Note that we do not 
assume yet that consecutive samples from each branch are independent. In fact, for the SC/OR 
analysis we will not assume any kind of fading distribution nor correlation model, allowing their 
arbitrariness. For the SSC, analysis in the IID branches scenario we will allow for arbitrary 
fading distribution with any correlation model, however we will have to take the additional 
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assumption of independent samples over time for characterizing the specific case of SSC with 
InID branches. This specific no correlation over time behavior can be found in time-division 
multiple-access (TDMA) systems, where the time interval between samples, T T dma = L x T s 
is assumed to be significantly greater than the channel coherence time, where L is the number of 
TDMA communications [26]. The SNR of each time interval is to be evaluated at its beginning 
so the chosen combining strategy can be performed based on this information and the availability 
of the diversity sources. 

A. Selection Combining 

SC combiner chooses each time the branch with the highest SNR. This technique benefits 
from fitting for both coherent and noncoherent modulation schemes, since the output of the SC 
combiner is equal to the signal on only one of the branches, hence the coherent sum of the 
individual branch signals is not required [9, Section 9.8]. 

1) Colocated Selection Combining: As the diversity sources are provided by several antennas 
at the destination, these are all permanently available at the receiver side. Thus, at the beginning 
of each time interval the SNR on each branch is evaluated and the one with the greatest one is 
selected to demodulate the message. In terms of the output SNR, this scheme can be expressed 
as 


Z[n] = max(Zj[n]) (1) 

2) Opportunistic Relaying: This term is usually found in literature for the SC distributed 
version. In this scheme, the destination has only one antenna and achieves spatial diversity 
through several relays. All relays listen to the source and attempt to decode the information, 
but only the one with the highest overall SNR z i[n] retransmits to the destination. In [10], it 
is described a procedure which achieves to make retransmit to the destination only to the relay 
with the highest overall SNR in each time interval without loss of spectral efficiency. Thus, 
given the availability of all the diversity sources all the time at the receiver, the output SNR can 
be expressed as in (1) too. 
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B. Switch & Stay Combining 

SSC combiner chooses one of the diversity branches and stays with it as long as the SNR on 
it remains above a preset threshold T. When this occurs, the combiner switches to the another 
branch. As it does not use signals from different branches together at one time, this technique also 
benefits from fitting for both coherent and noncoherent modulation schemes as SC [9, Section 
9.9]. The distribution of the output SNR will be given in terms of the distribution of the SNR 
on each input branch, but it will also depend on the instant the switching is to be performed, 
after switching condition is met, and on the order in which the different branches are selected. 
Because of this, we will consider that the system switching order is defined by a circular list of 
eligible branches, and that two different switching modes will be taken into account, each one 
related to the colocated and distributed scenarios respectively: 

1) Colocated Switch & Stay Combining: This is the classic manner of understanding SSC, as 
it is regarded in [9, Section 9.9.1]. A instant switching (IS) policy, i.e. the switching is performed 
at the precise time interval when the switching event occurs, is appropriate for colocated diversity 
systems, such as a multiple antenna receiver, since the next diversity source to be switched to 
(i.e., the next receive antenna to select) will be available at that very interval. Attending to IS 
policy, the output SNR of a colocated SSC system ca be expressed as 

if Z s [n] < T then 
s — (s + 1) 

end if 

Z[n] = Z s [n\ 

where s G |0... A r — 1] represents the branch that the combiner is tracking; Z[n] is the 
combiner output SNR in the nth interval; T is the switching threshold; and (•) represents the 
modulo operation, modf ■, A r ), which performs the circular step forward through the list of N 
eligible branches. 

2) Distributed Switch & Stay Combining: We understand as a distributed SSC system as the 
one which achieves its diversity by being able to link through several relaying nodes as it is 
regarded in [11], [12]. A deferred switching (DS) policy, i.e. the switching is to be performed 
in the next time interval to the one on which the switching event actually occurs, is the suitable 
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manner to perform SSC in a distributed system, where only one branch is available on each 
time interval, i.e. destination only links through a relaying node in each time interval, and the 
switching must be requested after the switching event occurrence. That is the destination will 
not detect the drop in the selected branch SNR until it receives it and must wait for the next time 
interval to request the link through another relaying node. Attending to DS policy, the output 
SNR of a colocated SSC system ca be expressed as 

Z[n\ = Z s [n\ 

if Z s [n] < T then 
s — (s + 1) 

end if 

C. Dual-hop relaying in distributed scenarios 

In this subsection, we specialize our distributed configuration for a dual-hop threshold-based 
DF strategy. As in [8] or [10], each time interval Ts is divided in two slots: during the first slot, 
the relays are overhearing from the source; in the second slot, the source does not transmits and 
waits for the selected relay, based on the switched combining strategy decision, to retransmit to 
the destination. Specifically the investigated scenario operates as follows: for the first hop, when 
the received SNR at the relay is greater than a decoding threshold T DF , noise free decodification 
is assumed and a perfect reconstruction of the transmitted signal is sent from the relay to the 
receiver. Thus, the distribution of the SNR at the receiver side only depends on the second hop 
link. When the received SNR at the relay is below T DF , the relay does not decode and transmits 
nothing; thus the SNR at the receiver side will be identically 0. As stated before, we assume 
the noise and fading caused by the first and second hop to be independent. 

This mode of operation corresponds mathematically in terms of SNR to 

7 r 1 / 0 if Z iA n \ < T DF, 

Zi[n\ = < (2) 

^ Z l)2 [n] if Z t . i [n] > T FF . 

III. Derivation of LCR and AFD 

In order to analyze higher order statistics of the output SNR of switched diversity combining 
schemes, we will use the alternative approach used in [24]. Specifically, the LCR and the AFD 
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can be obtained as: 


N z {u) = 


Az(u) = 


F z (u ) - F z (u,u) 


Ts 


F z (u) 


( 3 ) 

( 4 ) 


N z {uY 

where N z (u ) is the LCR of the stationary process Z[n] through the level u; A z (u ) is the AFD 
of Z[n\; F z {-) is the CDF of Z\n]; F z (-, •) is the bivariate CDF of the vector of two consecutive 
samples of Z[n\, Z = [Z[n], Z[n + 1]]; and T s is the time interval between samples. 

As (3) and (4) show, only the univariate and bivariate CDF of the output are needed to calculate 
its LCR and AFD. From now on, we will focus on obtaining the analytical expressions for these 
CDFs. 


A. Selection Combining 

As stated in section II-A, the output SNR can be express in terms of the per-branch SNR the 
same way for the colocated SC scenario or the distributed OR scenario. Thus, the expressions 
for the univariate and bivariate CDF are the same. Specifically, we can write the SC/OR output 
SNR univariate CDF in terms of the per-branch SNR CDFs as 


N-l 

F z (u) = Pr{max(Zj[n]) < u} = Pr {Z 0 [n] < u, Zi[n\ <u... Z N [n} < u} = F Zi (u), (5) 

i =o 

where we have made use of the independence between fading links assumption and F z< (•) is 
the univariate CDF of the SNR on the ith branch. In the same fashion, we express the output 
SNR bivariate CDF as 


Fz(u{, u 2 ) = Pr{max(Zj[n]) < w 1 ,max(Zj[n+ 1]) < u 2 } 

= Pr{Z 0 [n] < Mi... Z N [n] < Ui, Z 0 [n + 1] < u 2 ■ ■ ■ Z N [n + 1] < u 2 } 

N-l 

= n f Zi {ui,u 2 ), 


( 6 ) 


i =0 


where F Zi (•, •) is the bivariate CDF of the SNR of two consecutive samples Z % [n ], Z, [n + 1] 
on the /1h branch. 


August 12, 2016 


DRAFT 



10 


Given the expressions in (5) and (6), the LCR and AFD of the output SNR of any SC/OR 
system can be easily computed when expressions of the per-branch SNR univariate and bivariate 
CDF are available. 


B. Switch & Stay Combining 


We will first calculate expressions for the CDFs of interest, in terms of the distributions of 
the SNRs on each branch. For the colocated and distributed SSC cases, we will analyze two 
different scenarios: 

. IlD-AC : IID fading with arbitrary distribution and arbitrary temporal correlation model. 

. InID-TI: InID fading with arbitrary distributions and temporal independence between its 
samples. 

In general, the univariate and bivariate CDF of the SNR level at the output of the combiner 
can be expressed as 

N—l 

Fz{u) = ^2 p i F z\s=i(u), (7) 

i =0 


N -1 

F z (u 1 ,u 2 ) = yP l F^ i {uuU 2 \ (8) 

i =0 

where Pi is the probability of the combiner being tracking the zth branch at any given time 
and Fz\ s =i(u), F z \ s=i (ui, u 2 ) are respectively the univariate and bivariate CDF of the output 
SNR level, conditioned on the combiner being tracking the ith branch at the beginning of the 
observation interval. 

For the IID-AC case, it is trivial that P, : = TV -1 , F z \ s =i{u ) = F z \ s =j{u ) and F z \ s =i(ui, u 2 ) = 
F z \ s =j(ui,u 2 ) Vi,j e [0...7V — 1], 


For the InID-TI case, an expression for P, can be found in [25, Equation 2] as 

( n -i x -i 

"Kw 


( 9 ) 


Unfortunately, to the best of our knowledge, there is not an expression for P, in a SSC InID 
scenario which allows for correlation over time on each branch. 

Now we proceed to find the expressions of F z \ s =i(u), F z \ s =i(ui, u 2 ) for each of the investigated 
switching modes. 
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1) Colocated SSC: Assuming IS is performed, the combiner output SNR univariate CDF, 
conditioned on the combiner being tracking the zth branch at the beginning of the observation 
interval, reduces to an expression similar to [9, Eq. 9.270], 


Fzi(T) x F z (u) 


for u < T, 


Fz\s=i{^) \ 


( 10 ) 


Fz,(T) * F Zlw> (u) + F Zl (u) - F Zi (T) for u > T. 

On the other hand, the combiner output SNR bivariate CDF expression is divided in four 
intervals, ui,u 2 < T; u 2 < T < u\, U\ < T < u 2 , and ui,u 2 > T. By defining 


^n(ui,u 2 ) = ( 


F Zi {T,u 2 ) x F z (ui,T) 


for N = 2, 


(ID 


_ F Zi (T) x F z<i+1) (u u T) x F Z(i+2) (u 2 ) for TV > 3, 

the output SNR bivariate CDF corresponds to 




for ui,u 2 < T, 


^n(ui,u 2 ) 

for u 2 < T < ui, 

+ F Z{ . +1) (u 2 ) x \^F z .(ui,T) — F z .(T,T)j 


F z \ s =i(ui,u 2 ) 


^n(ui,u 2 ) 

+ F z . (T) x ^F Z{i+1) (u 1 ,u 2 ) — F Z{i+1) (ui,T)^ 


for u\ < T < u 2 . 


( 12 ) 


^n(ui, u 2 ) + F Zi (ui,u 2 ) - F Zi (ui,T) 

- F Zi (T, u 2 ) + F z .(T, T) 

, s for u\,u 2 >T. 

+ F Z{i+1) (u 2 ) x \ F z .{ui,T) — F Zi (T, T)j 

+ F z .(T) x ^F Z(i+1) (ui,u 2 ) - F Z{i+1) (ui,T)j 

We observe in (11), (12) that the bivariate CDF is different for N = 2 and N > 3; hence, 
the number of branches must be regarded when calculating statistics whose observation interval 
allows for the occurrence of two switchings. The expression must differ between a two-branch 
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combiner, N = 2, and a combiner with three or more branches, N > 3, since in the two-branch 
combiner case, after the second switch, the combiner returns to the original branch. 

2) Distributed SSC: Assuming DS is performed, the combiner output SNR univariate CDF, 
conditioned on the combiner being tracking the zth branch at the beginning of the observation 
interval, can be expressed as 

Fz\s=i(u) = F Zi (u). (13) 

The output CDF turns out to be the same CDF of the selected branch SNR, as the combiner 
outputs this branch regardless of its SNR level. The SNR level of the selected branch only 
determines the selected branch on the next interval. 

In this case, the bivariate CDF only has two regions and does not depend on the number of 
branches, since only one switch can occur during the observation interval of this second order 
statistic: 


F zX u i) x F z {i+1) (u 2 ) 


F z\s=i(u 1 ,u 2 ) = < 


for ui < T, 


F Zi(ui,u 2 ) — Fzi(T,u 2 ) + F Zi{T) x F Z(i+1) (u 2 ) for ui>T. 


(14) 


C. Distributions on each branch 

We now settle the expressions of the univariate and bivariate CDF of the SNR on each branch, 

FzA') and F Zi(-,-)- 

In the colocated diversity scenarios, this reduces to CDFs of each branch SNR. In the dis¬ 
tributed scenarios, we assume a DF relay network operating as indicated in section II-C. 

Considering that the SNR distribution on each hop is known, the overall SNR CDF on each 
branch can be expressed as 


Fz,M =Fz u (T DF ) + F Z( » x (1 - F z JT df )), (15) 

where F Zi , ( •) and F Zi 2 (•) are the CDFs of the SNR on the first and second hop of the z'th 
branch, respectively. 
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The bivariate CDF of each branch in the distributed scenario can be expressed as 

F Zi {u\,u 2 ) =F Zi1 (T df ,T df ) + {F Zi2 {ui) + F Zi2 {u 2 )) x [F Zil {T DF ) - F Zil (T DF ,T DF )) 
+ F Zi2 (ui, u 2 ) x (l + F Zi l (T DF ,T DF ) — 2 x F Zil (T DF )^, 


(16) 


where F Zil (-, •) and F z 2 (•, •) are the bivariate CDF of the SNR on the first and second hop of 
the ith branch, respectively. 

We must recall that when assuming independence between samples of a branch, every F z .{u \, u 2 ) 

F Zi (u i) x F Zi (u 2 ) 


IV. High Mean SNR Asymptotic Behavior 

While the approach discussed in section III is general and exact, it does not offer a solution 
when the bivariate CDF of the fadings involved in the scenario is not available. For this reason, 
and to provide some insight about the impact of the scenario parameters on studied higher order 
statistics, we proceed to develop simpler expressions which, only requiring the knowledge of 
the univariate PDF, describe considerably accurately the behavior of these higher order statistics 
in scenarios where the mean SNR is much greater than the level u. 

A. LCR in High Mean SNR Scenarios 

Let us note that the univariate and bivariate CDF of any fading SNR in addition to u is 
also composed of a parameter which we had omitted before for the sake of compactness: the 
mean SNR, O. Moreover, the correlation coefficient, p, has also been omitted from the bivariate 
CDF. Thus, F z (u ) = F z {u\Vt) and F z (ui,u 2 ) = F z (ui, u 2 \fi,, p). We define the normalized 
SNR, Z[n\ = Z[n]/Q, such as the expectation of Z is 1, and its CDFs F z (u), F z (ui, u 2 \p). By 
defining the normalized u, u = u/Vl we can rewrite the LCR expression in (3) as 

n z («) = (17) 

1 S 

Hence, if equation (17) is infinitely differentiable in u, the LCR for high mean SNR environ¬ 
ments can be expanded as a Taylor series at u — 0, 

1 ^u n d n . , i 

^(«) = - F z(«,M|p))| fi=0 - ( 18 ) 

b n =0 
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Typically, the univariate and bivariate CDF of a positive random variable is null at 0. Thus, 
the first non-null term of the LCR Taylor series expansion at u = 0 is enough to characterize 
the behavior of the LCR when asymptotically u /Q —> 0. Nevertheless, we have introduced 
distributed scenarios with DF relaying where the end-to-end SNR can be identically null with 
a finite probability, making theirs univariate and bivariate CDF non-null at 0. In this case, the 
asymptotic value of the LCR for low values of u is immediately the difference between the 
univariate and bivariate CDF divided by T s . In the rest of cases, it is clear that the LCR tends 
to be null for low values of u, but we attempt to know how it approaches to 0. Thus, we study 
now the behavior of the nth derivatives of the univariate and bivariate CDF 

Using equation [27, (0.410)], the nth derivative of any univariate CDF can be expressed as 


rjn J(n- 1 ) 

= *<^ (s) ' <19) 
where fz(u) is the PDF of the normalized fading SNR, whereas the nth derivative of the 

bivariate CDF can be expressed as 


d n ( r U - ^(n-l) " ^ ^k \ 

— F 2 {u,u\p) = 2( / fz 2 \zAMZi = u)dz 2 x —— iy /^(M) + ^$ n , fc (n,p)— f 2 (u) I, 

v>/0 fc =0 ' ' 

( 20 ) 

where f = ri) is the PDF of a sample of the normalized SNR conditioned to the 
value of the previous sample, and ^ n k (u,p) gathers together functions which multiply lower 
order derivatives of the PDF to compose the expression of the nth derivative of the bivariate 
CDF 

Note now that, if we assume the continuity of conditioned PDF around 0, the integral expres¬ 
sion multiplying the highest order derivative of the PDF involved in (20) for u = 0, becomes 0 
due to the limits of integration. Hence, in the inspection process in search of the first non-null 
term, we will evaluate at u = 0 the successive derivatives of the PDF until we locate the first 
non-null. At this point, we will have the value of the univariate CDF term, whereas we will 
have ensure that the bivariate CDF term is still null, since the first non-zero derivative term is 
multiplied by zero and the lower order derivative terms have been proved to be zero. 

If we call m to order of the first non-null derivative of the LCR at u — 0, we can express the 
LCR in high mean SNR environments as 
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Nz(u) |u_>o 


1 S m -Vf 2 {u) 


m\Tc du( m ~ 1 '> 


u 


( 21 ) 


w=0 


which we normali z e to the sampling period, Nz(u ) = Nz(u)Tg, and express in logarithmic 
scale to find that 


log(JV z (u))l«-+o ~ m log(fi) + lQ g(^y rf ' (22) 

the LCR in a high mean SNR scenario tends to a slope-intercept form where m, order of the 
first non-null derivative, determines the slope of the straight line asymptote; whereas the value 
of this mth derivative shifts left the asymptote. Thus we have obtained an simple expression for 
the asymptotic behavior of the LCR of a sampled process where only the knowledge of the PDF 
is required. 

B. AFD in High Mean SNR Scenarios 

After the asymptotic LCR analysis, recalling (4), assuming once again that the univariate and 
bivariate CDF of the process is null at 0, we can express the asymptotic behavior of the AFD 
as 


Az{u)\u-+o ~ Ts, (23) 

since the LCR asymptotic behavior in this case is the same to the univariate CDF asymptotic 
behavior. From this expression we observe that, for high mean SNR environments, the AFD 
tends to last the duration of a single sample of our discrete process, i.e. such fadings are so 
improbable that, in case of occurrence, its mean duration is the minimum possible. 

Whereas, when the univariate or bivariate CDF are non-null at 0, this A z (u) | u /n_». 0 is higher 
with direct proportion to how much probability is condensed at 0. For that case, we can write 

A zM\™ ~ Ts Fz (0) - F z (0,0)‘ (24) 
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V. Numerical Results 


In this section, we use the derived expressions for the higher order statistics to evaluate 
their behavior in different scenarios: SC, OR, colocated SSC and distributed SSC in Rayleigh, 
Nakagami-m and Hoyt fading environments. For the sake of simplicity, we plot the exact curves 
for Rayleigh fading scenarios, whose bivariate CDF follows a simple expression in terms of 
the Marcum Q-function [9, Eq. 6.5], verify them with simulation points, and compare them to 
the asymptotic straight line expression. While for the Nakagami-m and Hoyt scenarios we plot 
the low u /Q asymptotic straight lines and compare them with simulated points. Note that exact 
calculation of other fading distributions can be easily considered by plugging readily available 
expressions for the bivariate CDFs. We assume that the underlying Gaussians related to each 
probability distribution model of each fading link experience the temporal correlation model 
proposed by Clarke [28, Chapter 1]. Thus, the SNR correlation of each link is defined by the 
coefficient 


P= co ^Mt±M = \ M 2, fDTs )\\ 


(25) 


var (Z) 

where cov(-, •) is the covariance, var(-) the variance, J 0 (-) is the Bessel function of the first kind 
and order 0 and f F is the Doppler frequency of the channel. For the following results we use 
the values f F Ts = 0.05, 0.2 and 0.38, which lead to p « 0.95, 0.41 and 0. 

We use the levels u, T and T DF in a normalized manner, u = u/Q, T = T/Q, T DF = T DF /Q, 
where H is the mean SNR on each of the considered IID fading links. 

In Fig. 2, we show the normalized LCR of the output SNR in colocated diversity SC combiners 
with different numbers of antennas affected by IID Rayleigh fading. It is shown how high 
correlation between samples of our discrete processes leads to lower LCR, as it could have been 
forecast. We also observe how exact curves tend to approach to the asymptotic straight line for 
low u and that, the lower is the correlation, the higher values of u for which the asymptotic 
straight line accurately approximates the exact curve. In this scenario the increase of number of 
colocated antennas N makes the asymptotic behavior be steeper. Finally, the markers on these 
curves, and on all the following curves, correspond to simulations points which validate our 
analysis. 

Fig. 3 shows simulation results of normalized LCR of the output SNR in 2- and 4-antenna SC 
receivers, but IID Hoyt fading is assumed this time. We plot simulated curves to show that our 
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u( dB) 


Fig. 2. Normalized LCR of the output SNR in colocated SC combiners with 2 and 4 antennas in IID Rayleigh scenario. 


closed-form straight line asymptotic analysis is accurate for low u. We observe that the value of 
q, the shape parameter of the Hoyt distribution, moves the asymptote location, but not its slope. 
An explanation for this can be found if we take the Hoyt normalized power PDF, then we attend 
to (5), and apply the asymptotic analysis described in previous section. We obtain that 
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Fig. 3. Comparison of simulation results and asymptotic behaviour of normalized LCR of the output SNR in a colocated SC 
combiners with 2 and 4 antennas in q = 0.1 and q = 0.9 IID Hoyt scenario. 


log(iVz(u)) |u -> 0 ~ N log ( 1 + q ) + N log(w). (26) 

So the slope is given by the number of branches and the location by q and the number of 
branches as well. 

By setting q = 1, we obtain the asymptotes in Fig. 2, since the Rayleigh distribution equals 
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the q = 1 Hoyt distribution. 



Fig. 4. Normalized LCR of the output SNR in a colocated 3-antenna SC combiner in InID Rayleigh scenario. The mean SNRs 
on each InID Rayleigh are ft, fl/2 and S2/4 respectively 


We show a InID Rayleigh colocated SC scenario in Fig. 4. We impose to the first branch the 
mean SNR Q, which we also use to normali z e the level u, u — u/il, and then, fractions of it 
to the others, concretely 1/2 and 1/4. We can observe how the LCR is lower when the most 
powerful branch, i.e. the most selected, is also highly correlated in time. The expression for this 
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InID Rayleigh SC asymptote reduces to 


log(iV z (u))| u— >-0 


Nlog(u) — log ( XI 


N—l 


k =0 


(27) 


where (3 k is the ratio of the mean SNR on each branch and the mean SNR used to normalize 

u, /3k = Z k /Q, i.e. in this case f3 0 — 1, /3\ — 1/2, /3 2 = 1/4. 



Fig. 5. Normalized LCR of the output SNR in colocated SSC combiners in IID-AC Rayleigh scenarios and an additional 
InID-Tl scenario. 
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Fig. 6. Comparison of simulation results and asymptotic behaviour of normalized LCR of the output SNR in colocated SSC 
combiners with 2 and 3 or more antennas and T = —5 dB in m = 2 and m = 3 IID Nakagami-m with p = 0.95 scenarios. 


We plot on Fig. 5 the normalized LCR of the output SNR of colocated SSC receivers with 
N = 2 and N > 3 antennas in IID Rayleigh fading conditions with switching thresholds T — —5 
dB and 0 dB. There are exact curves for different correlation cases, whose validity is one more 
time verified by simulation points. Their straight line asymptotic behavior for low u is exhibited. 
The different number of antennas in this kind of combiners does not change their asymptote, 
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Fig. 7. Outage probability and Average Outage Duration (AOD) (outage level considered is 10 dB) against the chosen switching 
threshold T in colocated SSC combiners with N = 2 oi N > 3 IID Rayleigh diversity sources with per-branch mean SNR 
= 10, 20 dB and p = 0.95. 


while the switching threshold T changes its location but not its slope. The IV-dependence is 
more noticeable in the central parts of the graph near T. This is because, as we have observed 
in the asymptotic analysis, only the univariate CDF takes part for low u and this one is the 
same no matter N, see (10). And for high u, attending to (12) common new terms add up to 
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the N-dependent term defined in (11). Moreover, we made a difference on the bivariate CDF 
regarding N because the correlation effect this CDF experiences when returning to the original 
diversity branch after a second switch or switching to a third branch, see (11). That is why we 
cannot see the N = 2 and p — 0 curve, the N > 3 and p = 0 curve is the same and has been 
plotted over the N = 2 one since there is no difference when different N is available at the 
receiver and there is no correlation over time on each fading branch. In addition we also plot 
on Fig. 5 a SSC InID-TI example, concretely the dual branch SSC where on one branch we 
have Rayleigh fading and Hoyt (q = 0.3) on the other, and we verify its validity with simulation 
points as well. 

In Fig. 6 the results of colocated SSC receivers in IID Nakagami-m environment are plotted. 
We see how, this time, the m parameter affects on the asymptote slope and the switching threshold 
to the location. Nevertheless, the number of diversity sources does not change the asymptotic 
behavior for this kind of combination. An explanation for this can be found if we take the 
Nakagami-m normalized power PDF, then we attend to (10), and apply the asymptotic analysis 
described in previous section. We obtain that 

log(N z (u))\ u ^ 0 ~ log(——— —-^F z (T\m,tt)) +m\og(u). (28) 

VI (m) m\ J 

So the slope is given by the parameter m and the location by the switching threshold and m 
as well. 

By setting m = 1, we obtain the asymptotes in Fig. 5, since the Rayleigh distribution equals 
the m — 1 Nakagami-m distribution. 

Fig. 7 goes deeper on the performance analysis of SSC systems by adding the study of second 
order statistics of the outage. We see how the choice of a proper switching threshold in SSC 
systems is important to reduce the probability of an outage event to occur (left scale on the 
graph), and that when the threshold equals the outage level the outage probability minimizes. 
In addition, we observe the Average Outage Duration (AOD) (right scale on the graph) and 
realize that its minimum is not necessarily reached for the same threshold which minimizes 
the outage probability. This fact is important to consider for applications where latency of the 
communication is critical, such as real time communications, where the outage duration bounds 
the latency limits. Finally, the most important behavior we observe from this curves is that the 
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average outage duration is higher for the dual branch SSC combiner than for the SSC combiner 
with a third branch or more, which refutes the popular believe that SSC systems do not benefit 
from more than two diversity branches [25]. 


x 


£ 



Fig. 8. Normalized LCR of the output SNR in IID Rayleigh dual hop dual branch (N = 2) Opportunistic Relaying scenario. 


Fig. 8 shows the normalized LCR of the SNR at a single antenna receiver from a dual hop 
dual branch threshold based DF OR scheme when all the links involved in the configuration 
undergoes IID Rayleigh fading. We plot exact curves verified by simulation points. For this 
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u(dB) 


Fig. 9. Normalized AFD of the output SNR in IID Rayleigh dual hop dual branch (N = 2) Opportunistic Relaying scenario. 


configuration the asymptotic behavior is a horizontal straight line due to threshold based DF 
mode of operation, which condenses a finite amount of probability at Z = 0. Attending to (5), 
(6), (15) and (16), in the IID case, we can write the expression for the horizontal asymptote, 


N z (u )\^o ~ (Fz(T df )) N - (Fz(T df ,T df \ P )) N , 


(29) 
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where F z {z) and F z (z, z\p) represent the univariate and bivariate CDF of the IID fading on 
each branch. 

We observe that in this particular scenario the correlation affects the asymptotic behavior, in 
particular on its location: low correlation rises the asymptote. Another parameter which changes 
the asymptote location is the decoding threshold Top. 

In Fig. 9 we plot the normalized AFD for the same scenarios as Fig. 8. The AFD does not 
tend to 1 for low u as we had foreseen in previous section due to the condensed probability at 
Z = 0 caused by the threshold based DF relays. Higher Tof means higher probability of the 
relays not forwarding the signal, so the AFD rises. The higher temporal correlation also makes 
the AFD rise. These two phenomenona are evidenced in the figure. 

In Fig. 10, it is shown the normalized LCR of the SNR at a single antenna receiver from a 
distributed SSC dual hop threshold based DF relaying scheme when all the links involved in 
the configuration undergoes IID Rayleigh fading for any N. We plot exact curves verified by 
simulation points. For this configuration the asymptotic behavior is again a horizontal straight 
line due to threshold based DF mode of operation. Concretely, attending to (13), (14) and (15), 
in the IID case, 


Nz(u)\u->o ~ F z {T df ){1 — F z (Tdf )). (30) 

We observe here that the location of the horizontal asymptotes only depends on the decoding 
threshold Tof, but not on the per-branch correlation. In fact, this scheme exhibits little difference 
for different correlations. 


VI. Conclusion 

We derived novel expressions for the LCR and AFD of communication systems based on 
the concept of switched diversity diversity to achieve full-order spatial diversity. Threshold 
based techniques lacked from this analytical characterization for these statistics because of the 
inherent discontinuity of the random processes of interest; we circumvented the limitations of 
Rice’s approach by using an alternative formulation for sampled random process. We achieved to 
characterize SC and SSC techniques, which are now again subject of interest since the introduc¬ 
tion of cooperative diversity, in classic colocated diversity and already mentioned cooperative 
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Fig. 10. Normalized LCR of the output SNR in distributed SSC combiners in IID Rayleigh scenario (p = 0.05). 


distributed diversity scenarios in a discrete time fashion, which matches the actual inherent 
implementation of these techniques. Furthermore, the SSC schemes lacked of LCR and AFD 
analytical characterization until now to the best of our knowledge. We also provided another 
simple manner to elude the discontinuity when analyzing the higher order statistics of threshold- 
based DF relaying systems. 

Our analysis holds for an arbitrary fading distribution either when assuming colocated or 
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distributed diversity branches. In the process of computing the LCR and the AFD, we obtained 
analytical expressions for the univariate and bivariate CDFs of the output SNR. We also studied 
the effect of temporal correlation in these scenarios. Moreover, we introduced a common analysis 
for the asymptotic LCR and AFD in high mean SNR environments which allows to characterize 
the behavior of the higher order statistics only requiring the availability of the univariate PDF 
which describes the fading random process. We also discussed the implications of using a number 
of branches in SSC systems larger than the conventional recommendation of two for wireless 
systems finding interesting results about how more than two branches improves the average 
outage duration. Finding that a third branch in SSC systems reduces the average outage duration 
of these. 
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